Derivation and Implementation of the Gradient of the Râ•fi7 Dispersion Interaction in the Effective Fragment Potential Method

The dispersion interaction energy may be expressed as a sum over R−n terms, with n ≥ 6. Most implementations of the dispersion interaction in model potentials are terminated at n = 6. Those implementations that do include higher order contributions commonly only include even power terms, despite the fact that odd power terms can be important. Because the effective fragment potential (EFP) method contains no empirically fitted parameters, the EFP method provides a useful vehicle for examining the importance of the leading R−7 odd power term in the dispersion expansion. To fully evaluate the importance of the R−7 contribution to the dispersion energy, it is important to have analytic energy first derivatives for all terms. In the present work, the gradients of the term E7 ∼ R−7 are derived analytically, implemented in the GAMESS software package, and evaluated relative to other terms in the dispersion expansion and relative to the total EFP interaction energy. Periodic boundary conditions in the minimum image convention are also implemented. A more accurate dispersion energy contribution can now be obtained during molecular dynamics simulations. ■ INTRODUCTION Dispersion derives from the intermolecular forces that arise from the interaction of induced multipoles between atoms and molecules. The dispersion interaction between two molecular species can be expressed as = + + = + + + E E E E C R C R C R ... / / / ... disp 6 7 8